Resonant quantum kicked rotor with two internal levels 
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We develop a system consisting of a quantum kicked rotor with an additional degree of freedom. 
This models a single two-level atom with internal ground and excited states, and it is characterized 
by its quantum resonances with ballistic spreading and by the entanglement between the internal 
and momentum degrees of freedom. These behaviors establish an equivalence between our model 
and the usual quantum walk on the line. 
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I. INTRODUCTION 

Advances in technology during the last decades have 
made it possible to obtain samples of atoms at tem- 
peratures in the nK range pj (optical molasses) using 
resonant or quasiresonant exchanges of momentum and 
energy between atoms and laser light. The experimen- 
tal progress that has allowed to construct and preserve 
quantum states has also opened the possibility of build- 
ing quantum computing devices [^-Q and has led the 
scientific community to think that quantum computers 
could be a reality in the near future. This progress has 
been accompanied with the development of the interdis- 
ciplinary fields of quantum computation and quantum 
information. In this scientific framework, the study of 
simple quantum systems such as the quantum kicked ro- 
tor (QKR) d, Q and the quantum walk (QW) [H may 
be useful to understand the quantum behavior of atoms 
in optical molasses. 

The QKR is considered as the paradigm of periodically 
driven systems in the study of chaos at the quantum level 
. This system shows behaviors without classical equiv- 
alent, such as quantum resonance and dynamical local- 
ization, which have posed interesting challenges both in 
theoretical and experimental @ terms. The occurrence 
of quantum resonance or dynamical localization depends 
on whether the period of the kick T is a rational or irra- 
tional multiple of Att. For rational multiples, the behavior 
of the system is resonant while for irrational multiples the 
average energy of the system grows in a diffusive manner 
for a short time and then the diffusion stops and localiza- 
tion appears. From a theoretical point of view the two 
types of values of T determine the spectral properties 
of the Hamiltonian. For irrational multiples the energy 
spectrum is purely discrete and for rational multiples it 
contains a continuous part. Both resonance and local- 
ization can be seen as interference phenomena, the first 
being a constructive interference effect and the second 
a destructive one. The QKR has been used as a theo- 



retical model for several experimental situations dealing 
with atomic traps 
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and is a matter of permanent 

attention pK26|. 

The quantum walk has been introduced [1, [27l-l33l as 
a natural generalization of the classical random walk in 
relation with quantum computation and quantum infor- 
mation processing. In both cases there is a walker and a 
coin; at every time step the coin is tossed and the walker 
moves depending on the toss output. In the classical 
random walk the walker moves to the right or to the left, 
while in the QW coherent superpositions right/left and 
head/tail happen. This feature endows the QW with out- 
standing properties, such as the linear growth with time 
of the standard deviation of the position of an initially 
localized walker, as compared with its classical counter- 
part, where this growth goes as i^/^. This has strong im- 
plications in terms of the realization of algorithms based 
on QWs and is one of the reasons why they have re- 
ceived so much attention. It has been suggested [s^ that 
the QW can be used for universal quantum computation. 
Some possible experimental implementations of the QW 
have been proposed by a number of authors 0, |36| - 
Soll. In particular the development of techniques to trap 
samples of atoms using resonant exchanges of momentum 
and energy between atoms and laser light may also pro- 
vide a realistic frame to implement quantum computers 
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A parallelism between the behavior of the QKR and 
a generalized form of the QW was developed in Refs. 
[21, 22] showing that these models have similar dynam- 
ics. In those papers, the modified QW was mapped into a 
one-dimensional Anderson model [42|], as had been pre- 
viously done for the QKR [43|. In the present paper, 
following the work of Saunders et al. [3,14^ we propose 
a modification of the QKR. We study some properties 
of this new version of the QKR and establish a novel 
equivalence between this new QKR and the QW. Essen- 
tially, the new QKR has an additional degree of freedom 
which describes the internal ground and excited states of 
a two-level atom. We call this new system the two-level 
quantum kicked rotor (2L-QKR). In this system the in- 
ternal atomic levels are coupled with the momentum of 
the particle. This coupling produces an entanglement be- 
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tween the internal degrees of freedom and the momentum 
of the system. 

The rest of the paper is organized as fohows, in the 
next section we present the 2L-QKR system. In the third 
section we obtain the time evolution of the moments. In 
the fourth section the entanglement between the internal 
degrees of freedom and momentum is studied. In the last 
section some conclusions are drawn. 



II. TWO-LEVEL QUANTUM KICKED ROTOR 

We consider a Hamiltonian that describes a single two- 
level atom of mass M with center-of-mass momentum 
described by the operator P. Its internal ground state 
is denoted by the vector \g) and its excited state by the 
vector |e). The internal atomic levels are coupled by two 
equal-frequency laser traveling waves with a controllable 
phase difference. Following |4J| , after a shift of the energy 
values, the 2L-QKR Hamiltonian can be written as 



pi 



H = — + h^\e){e\ 



-K6T{t)cos{kLz){\e){9\ + \g){e\). 



(1) 



Here A is the detuning between the laser frequency and 
atomic transition frequency. K is proportional to the 
Rabbi frequency and we shall refer to it as the strength 
parameter. 



drd) ^^5{t~~ nT) 



(2) 



n=0 



is a series of periodic Dirac's delta applied at times t — 
nT with n integer and T the kick period, z is the operator 
of the atom's center of mass position. Finally, /cl is the 
laser wave-vector magnitude along the z direction. 

Unlike the QKR, in the 2L-QKR the conjugate position 
and momentum operators have discrete and continuous 
components, i.e. 



It must be noted that Eq.([5]) does not depend on the 
operator I and therefore /3 is a preserved quantity. Then 
if the initial condition belongs to a subspace correspond- 
ing to a well defined eigenvalue of the dynamics is 
such that the system remains in said subspace and the 
evolution of the system will be only determined by the 
conjugate operators 9 and k. Therefore we may restrict 
ourselves to the study of the evolution constrained to 
a subspace corresponding to a given eigenvalue of /3. In 
this case the composite Hilbert space for the Hamiltonian 
Eq.(l5]) is the tensor product 'Hs ®'Hc- Tis is the Hilbert 
space associated to the discrete momentum on the line 
and it is spanned by the set {|fc)}. T-Lc is the chirality (or 
coin) Hilbert space spanned by two orthogonal vectors 
|e)}. In this composite space the system evolves, at 
discrete time steps t G N, along a one-dimensional lattice 
of sites fc G Z. The direction of motion depends on the 
state of the chirality. Taking this into account it is clear 
that the Hilbert space of the 2L-QKR (with the preced- 
ing restriction) is identical to that of the usual QW on 
the line. 

The evolution of the system is governed by the Hamil- 
tonian given by Eq.([5|), so that, as is the case for the 
usual QKR, the unitary time evolution operator for one 
temporal period T can be written as the application of 
two operators, one representing the unitary operator due 
to the kick and another being the unitary operator of the 
free evolution [ill 

where ax is the Pauli matrix in the x direction. 
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(27r/ -f 9) 



P^hkL{k + l3) 



(3) 



(4) 



where the eigenvalues of I and k are integers and the 
eigenvalues of 6* G [— tt, tt) and the eigenvalues of the 
quasimomentum j3 £ [—1/2,1/2). It is important to 
point out that the operator /3 commutes with both k 
and 9. Using Eqs.(|3l4l) to substitute z and P in Eq.^ 
yields 



H 



hkL{k + !3) 



ftA|e)(e| 



2M 

+K5T{t)cos{9){\e){g\ + \g){e\). (5) 



The unit operator Eq.Q in the momentum representa- 
tion and in the chirality base {|e), \g)} has the following 
shape 



<5fc-j 21 + 1 Sk^j 21 



(9) 



e Ok-j 21 e Ok-j 21+1 



where 



6kj is the Kronecker delta, I is an integer number and 
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The wave- vector in the momentum representation can be 
expressed as the spinor 



(12) 



A;— — oo 2 ^ ^ 

where /3 is the value of /3' for the chosen subspace and 



are the upper and lower components that correspond to 
the left and right chirality of the QW. 

The discrete quantum map is obtained using Eas. (|9ll2p 



j = -oo 



bk+pit + T) 



(14) 



The dynamical evolution of the system up to i = nT is 
obtained applying the above rule Ea. (fT4| n times. 



A. Resonance r = 27r in the /? = subspace with 

A = 2rmv 

In this subsection we solve analytically the evolution of 
the system given by the map Ea. (|14p . We consider here 
the principal resonance r = 27r in the subspace (3 = 0. 
Due to the quasimomentum conservation the value of /3 
does not change. Therefore the accessible momentum 
spectrum is discrete and from now on the theoretical de- 
velopment is similar to that of the usual QKR in reso- 
nance. Additionally we choose A = 2m7r with m integer 
in order to obtain the wave function analytically. We 
will show afterwards, using numerical calculation, that 
the qualitative behavior will be similar for arbitrary A. 
With these conditions the matrix of Eq.® only depends 
on j — k. In order to simplify the notation we define 



\ Oka-kb 21 + 1 (>ka-k,, 21 

(15) 

Using Eq. ([T^ the initial condition is connected with the 
wave function at the time t = nT by the equation 



afe„ {nT) 
bkMT) 



EEE 

kn — l kn — 2 kn — 3 



- E E E ^fe-fe—i (.'^Wk^-ik^-2 ('«)•■• 

k2 ki ko 

■■■Uk.kA^)UkM^)(^'f^oy (16) 

where a^^ = UkoiO) and 6^,^ = 6fe„(0). 



Using the relation, 

^ C/fc„fc„_i('«l)C^fc„_ifc„_2(K2) = t/fc„fc„_2('«l + K2), 

(17) 

obtained in Appendix A, Eq. p^ is reduced to 



afc(nT) 

bkinT) J ^2^'^ •J,-kK<^<^)y,-k2iy^^,^ 



^ E*^ ^-^^-^ ('^'*) 1 



k 21+1 



i 



(18) 



where I is now an arbitrary integer number. 



B. Antiresonance t = 2n in the /3 = subspace with 

A = {2m + l)n 

We now find the time evolution of the wave function 
for A = {2m + l)7r. Eq.® shows that in this case the 
matrix U{/3 = 0)^^ satisfies the relation 

C^fc„fc„_i(K)f/fc„_ifc„_2(«^) = 4„ fc„_2-^: (19) 

where / is the identity matrix. This last expression to- 
gether with Eq. p^ imply that 



ak{nT) 
bk{nT) 



'n 21 + 1 



E^.(-)(^) 



(20) 



Then it is clear that the 2L-QKR shows a periodic be- 
havior when the parameters of the system take the values 
here considered. This behavior has no analog in the usual 
QKR since the parameter A does not exist in said sys- 
tem. Furthermore, it is interesting to point out that this 
anti-resonance occurs for r = 27r, value for which the 
usual QKR is in resonance and does not present periodic 
behavior. 



III. PROBABILITY DISTRIBUTION OF 
MOMENTUM 

The evolution of the variance, ct^ = 7712 — TOi, of the 
probability distribution of momentum is a distinctive fea- 
ture of the QKR in resonance. It is known that it in- 
creases quadratically in time in the quantum case, but 
only linearly in the classical case. In this section we study 
the evolution of the variance of the 2L-QKR, once again 
restricting ourselves to the /3 = subspace and taking 
T = 27r, which corresponds to the primary resonance of 
the usual QKR model. We will obtain the variance from 
the evolution of the first and second moments, defined 
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as mi{t) = ^i^kPk{t) and m2(t) = ^k^'^^k{t) respec- 
tively, where Pk{t) — \ak{t)\'^ + \bk{t)\'^ is the probabihty 
to find the particle with momentum p = hkLk at time t. 

We first consider the resonance defined by A = 2m7r. 
In this case we are able to calculate the first and second 
moments analytically using Eq. p8)) and the properties of 
the Bessel functions (see Appendix B), obtaining: 



TOi(n) 



m2{n) = 



E 



■mi(0), (21) 



(Kn)^ 



1+ y: ^{c^a^i2(o)+m 



+2J 



+Kn J2 (2j + m K^";i + a°,bf_,] 

j=-oc 

+m2(0), 



(22) 



where [a;] and S [a;] are respectively the real part and 
imaginary part of x. mi(0) and TO2(0) are the moments 
at time t = 0. These last equations show that the behav- 
ior of the variance = m2 — has a quadratic time 
dependence irrespective of the initial conditions taken. 

When A = (2n -|- l)7r, it was shown in the previous 
section that the 2L-QKR has a periodic dynamics and 
therefore the behavior of the statistical moments will be 
periodic as well. 

The case when A 7^ 2mr is cumbersome to solve analyt- 
ically, so we restrict ourselves to a numerical study. The 
evolution of the second_statistical moment was obtained 
for different values of A through numerical iterations of 
the map given by Eg. ([Til) . It was found, for all the con- 
sidered values of A 7^ 2mr, that the long-time behavior 
of the second moment (and therefore of the variance) is 
quadratic after an initial transient. The duration of the 
initial transient depends on the initial conditions and the 
value of A. This features can be appreciated in Fig.([T]). 
The figure shows the time evolution of the second mo- 
ment for the initial conditions |'I'(0)) = \k = Q)\g)- It 
can be appreciated that the second moment approaches 
a quadratic behavior after an oscillatory transient. It was 
found that the nearer the parameter A is to (2?! -f- 1)tt, 
the more pronounced this oscillation is. 



IV. ENTANGLEMENT 

In the context of QWs several authors [53-l66| have 
been investigating the relationship between the asymp- 
totic coin-position entanglement and the initial condi- 
tions of the walk. In order to compare the model con- 
sidered in this paper with the QW, we investigate the 
asymptotic chirality-momentum entanglement in the 2L- 
QKR. The unitary evolution of the 2L-QKR generates 
entanglement between chirality and momentum degrees 
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FIG. 1: The dimensionless second moment for A — O.QTtt as 
a function of the dimensionless time. 



of freedom. This entanglement will be characterized 
[54i . i65i] by the von Neumann entropy of the reduced den- 
sity operator, called entropy of entanglement. The quan- 
tum analog of the Gibbs entropy is the von Neumann 
entropy 



Sn{p) = -tr(plogp). 



(23) 



where p — \'^{t)) (^(i)| is the density matrix of the quan- 
tum system. Owing to the unitary dynamics of the 2L- 
QKR, the system remains in a pure state, and this en- 
tropy vanishes. In spite of this chirality and momentum 
are entangled, and the entanglement can be quantified 
by the associated von Neumann entropy for the reduced 
density operator: 



S 



-tr(pc log2 Pc), 



(24) 



where pc —tTk{p) is the reduced density matrix that re- 
sults from taking the partial trace over the momentum 
space. The reduced density operator can be explicitly ob- 
tained using the wave function Eg. dT^ in the subspace 
/3 = and its normalization properties 



where 



Pg{n) Q{n) 
Q*{n) P,{n) I ' 



Pg{n)= £ \ak{nT)\\ 

00 

Pe{n)^ ^ MnT)\\ 

00 

= E ak{nT)bl{nT). 

j = -CX} 



(25) 



(26) 



(27) 



(28) 
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Pe{n) and Pg{n) may be interpreted as the time- 
dependent probabilities for the system to be in the ex- 
cited and the ground states respectively. In order to 
investigate the entanglement dependence on the initial 
conditions, we consider the localized case, that is the ini- 
tial state of the rotor is assumed to be sharply localized 
with vanishing momentum and arbitrary chirality, thus 

i rSn = f O^o, (29) 

\hk{Q) ) l^cxp«(^sm§y ™' ^ ' 

where 7 G [0, tt] and Lp £ [0, 27r] define a point on the 
unit three-dimensional Bloch sphere. Eq. ([T8l) takes the 
following form 

Substituting Eq.dSOl) into Eqs. (|26|27|28p and using the 
properties of the Bessel functions, we obtain: 

- i [l-f Jo(2nK)cos7], (31) 
Pe(ri) = Jo(2nK)cos7], (32) 



and the asymptotic value of the entropy, S{n) — 5*0, is 

So = - A+ log2 A+ - A_ log2 A- . (37) 

For the initial condition — 11/2 and/or 7 = tt on the 
Bloch sphere, Q ^ and both eigenvalues are A± = 1/2. 
In this case the asymptotic entanglement entropy Eg. ([57)1 
has its maximum value = 1 . Finally, for sharply local- 
ized initial conditions with zero momentum. Fig |2] shows 
the dependence of the asymptotic entanglement entropy 
on the parameters Lp and 7. 




0.2 0.4 0.6 0.8 1 



Q(n) 



sm7 



cos ip — i sin pj^ (2nK)] 



(33) 



The eigenvalues of the density operator pc, Eq. ((25| . as a 
function of Pg{n), Pe{n) and Q{n) is 



1 



l±yi-4(p<,(n)Pe(n)-|g(n)|2) 



, (34) 



and the reduced entropy as a function of these eigenvalues 
is 



S{n) = -A+ log2 A+ - A_ log2 A_ 



(35) 



Therefore the dependence of the entropy on the initial 
conditions is expressed through the angular parameters 
ip and 7. This means that, given certain initial condi- 
tions, the degree of entanglement of the chirality and 
momentum degrees of freedom is determined. 

It is seen from Eas. (|31l32l33p that the occupation prob- 
abilities and the coherence Q tend to a certain limit when 
n — > 00. In this limit Jo(2nK) — >■ and both of the occu- 
pation probabilities tend to 1/2, irrespective of the initial 
conditions. However, in the asymptotic regime, depen- 
dence on the initial conditions is still maintained by Q, 
and therefore by the entropy as well. Thus, in the asymp- 
totic regime we have 



A4 



A4 



1 



[1 ± cos ip sin 7] 



(36) 



FIG. 2: The dimensionless entanglement entropy as a func- 
tion of the dimensionless initial conditions, see Ea. (|29p . The 
grayscale (color online) corresponds to different values of the 
entropy between zero and one. 



V. CONCLUSION 

We developed a new QKR model with an additional 
degree of freedom, the 2L-QKR. This system exhibits 
quantum resonances with a ballistic spreading of the vari- 
ance of the momentum distribution, and entanglement 
between the internal and momentum degrees of freedom 
only depending on the initial conditions. These results 
were established analytically and numerically for differ- 
ent values of the parameter space of this system that 
correspond to the primary resonance of the usual QKR 
model. The above two behaviors also characterize the 
QW on the line and hence establish again an equiva- 
lence between the QW and the 2L-QKR. This suggests 
that experiments that are related to each of two mod- 
els should also carry some kind of physical equivalence 
between them. We have found also that, although our 
system exhibits characteristics similar to those found in 
the usual QKR model, there are still novel features, such 
as the existence of the anti-resonance described in section 
II B, which have no analogue in the simple QKR model. 
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These characteristics of the 2L-QKR render the system 
as an interesting candidate for further study within the 
framework of quantum computation. 

We acknowledge stimulating discussions with Victor 
Micenmacher, the support from PEDECIBA and ANIL 



Appendix A 

Starting from Eg. dTSl) the following expression is ob- 
tained 

J2 Uk^''^ (^)Ukok, {^) = J2 ( % eI 

ki ki 

where 

-i2A r c , ~iA 

Oui 21 0i,2 21' + e 

E2 = e-'^^S,, 21 (1 - S,, 21') + e-'^ 2V (1 - 21) , 
Ea = e"*'^ 21' (1 - Si,^ 21) + 21 (1 - 21') , 

Ea = e"''^ {l~S^^2l) i'^ - 21') + Ki 21K2 21', 

and with vi = ki — ka, V2 — k2 ~ ki. In the above equa- 
tions, three different type of sums are involved, which can 
be carried out using the properties of the Bessel functions 
(Ref.'eV^, p. 992, Eq. 8.530). 



where 



(Al) 



El = e-"^(5,, 21 21' + e~"^ (1 - 5,, 21) (1 " 5,, 21') , 



Jk2-ki (k) Jfci-fco ('«) — Ji_i2 ; 



(A2) 



Jk2-kt (k) Jki-ko (k) Ski- 



kn 2/ 



[Jf,2 (2k) +Sk2ko] , 



(A3) 



^ Jk2-ki (k) Jki-ko (k) Ski- 



-,Sfj.2 2{l + l') [Jl-12 (2k) +5k2ko] 



(A4) 



where H2 = k2 — kg. Substituting the above equations 
into Eq. (|A1[) and defining p = I + I' 



J2^kik2il^l)UkokAl^2) 



-lA 



z^V^, (2«) S^2 2p{l + e-'^y 
z'^V^, (2k) (1 + e-^^) (1 - V.2p) 
i'^V^, (2k) (1 + e^^) (1 - S^2 2p), 
i^-J^2 (2k)V2 2p (l + e^^), 
Sk2ko (e-^^ - 1) , 
Sk2ko fe'"^ - 



Fi 
F2 
F3 
F4 
Gi 
G2 



Appendix B 

The probability Pk{n) of finding the system with mo- 
mentum fc at a time t — nT is obtained using Eq. ()18p . 



Pfe(n) = |afe(n)|2 + |6fe(n)|2 = 

where 

Iji = [Jk-j{nn)Jk-i{nK) + Jk-j{-nK)Jk-i{-nK)] 

and a^. and are given by the initial conditions of the 
system. To calculate the moments mi(n) and m2{n) we 
need the following sums 

00 

4'' = E ^^fe-j(«)-/fc-/('^) 



and 



k'^Jk-j{n)Jk-i{n) 

fc— — 00 



(B2) 
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A ^3 






G2 J 



(A5) 



1 

= y('^'^' - 2{^l3+'2 + Sij-2)) 

+iK[^{5ij+i +Sij-i)+j{Si-i+i ~Sij-i)] +fSji. 

(B3) 

Using these expressions together with Eq. (|Bip and the 
definition of the moments we obtain the first and second 
moments Eqs. (j21l22l) . 
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